YAO Jun, HE Jiao, WU Yongfei, et al. Energy optimization of electricity wholesale market considering carbon emissions trading and green power certificate trading mechanism[J]. Electric Power, 2022, 55(8): 190-195.
引言
近年来,一些发达国家和地区针对电力行业制定了一系列碳减排组合机制[1]。其中,碳排放权交易(简称碳交易)和绿色电力证书交易(简称绿证交易)制度组合已经在美国、加拿大、英国、日本以及欧盟等国家和地区广泛应用,并对推动清洁低碳发挥了重要作用[2]。在碳交易和绿证交易制度双重约束下,如何合理确定最优碳配额系数和绿证比例组合,从而实现电力行业碳减排目标,确保区域内发电商利益,是亟待解决的问题。目前,国内外主要集中于单一碳交易或绿证交易制度下电力市场均衡问题研究。文献[3-7]构建了电力市场均衡模型,分析了发电商的策略性行为。文献[8-9]采用可计算一般均衡模型(computable general equilibrium,CGE)、双重差分模型(differences-in-differences,DID)、合成控制法(synthetic control method,SCM)等评估了碳交易政策。文献[10-11]运用电力市场均衡模型、系统动力学模型等分析绿证交易制度对电力市场的影响。文献[12]分析了碳交易机制与可再生能源支持政策相互作用对电力市场的影响。文献[13]考虑了可再生能源补贴政策对电力市场的影响,构建了发电商利润最大化为目标的决策模型。文献[14]分析了碳市场、绿证市场以及电力市场的相互作用关系,但仅定性说明了碳市场和绿色市场对发电企业发电决策的影响。文献[15]将碳交易和绿色证书交易机制引入到发电权交易模型中,构建了基于区域净利润最大化的发电权交易模型,但没有考虑不同电源之间的发电权交易。
(3)各发电商的发电成本为二次函数,即式中: qk为第 k 个发电商的实际发电量;αk、βk、γk分别为第 k 个发电商的成本参数; Ck为第 k 个发电商的发电成本。可再生能源发电具有强烈的间歇性和波动性,其需要投入比火力发电更多的运维和旋转备用成本。本文假设不同技术类型常规火力发电机组的发电边际成本始终要低于可再生能源发电边际成本[20]。(4)电力批发市场中电力供需始终处于平衡状态,电力批发价格与电力总需求量之间满足线性函数关系。电力批发价格pe的逆需求函数可以表示为式中:A、B分别为批发价格系数,且均大于0。(5)碳交易价格由市场供求关系决定,满足的线性函数关系为[21]式中:为碳交易价格;a 、 b 均为大于0的常量;eh,i为第i个火力发电商的实际碳排放量;为管制机构确定的区域范围内免费碳配额上限。同样地,假设绿证交易价格由市场供求关系决定,并且满足如下函数关系式[22-23]。式中:ρ为绿证交易价格; ξ 、τ分别为大于零的常量;λ为各发电商须义务完成的可再生能源发电量比例。由于第i个火力发电商发出的是非可再生能源电量,故该发电商需购买数量为 λ 的绿证以满足政府配额考核要求。第j个可再生能源发电商将数量为 (1−λ)qr,j的绿证在绿证市场上售出。2.2 模型建立在综合考虑碳交易和绿证交易制度后,由于火力发电商和可再生能源发电商发电成本、碳交易和绿证交易成本上存在差异,故不同类型发电商的净收益函数并不相同。火力发电商的净收益主要包括发电净收益、买入碳配额支付的成本或卖出碳配额获得的收益、支持绿色电力发展的费用。第i个火力发电商自身净收益函数 Rh,i为式中:为第i个火力发电商实际买入或卖出碳配额的数量;Ch,i为第i个火力发电商的发电成本。可再生能源发电商除了在电力市场中获得一部分发电收益外,还能在绿证市场中售出一部分绿证获得收益。第j个可再生能源发电商的净收益 Rr,j为式中:为第j个可再生能源发电商实际卖出的绿证数量;Cr,j为第j个可再生能源发电商的发电成本。从发电资源合理配置的角度,需要实现区域内各发电商总收益最大化来确定发电商的上网电量[24-25]。区域内各发电商总收益R为在考虑碳交易和绿证交易制度双重约束下,各发电商总收益最大化函数受到5个约束条件,具体如下。(1)火力发电商实际买入或卖出的碳配额量不超过自身需要买入或能够卖出的碳配额量;(2)各火力发电商实际碳排放总量不超过管制机构规定的免费碳配额上限;(3)各火力发电商实际买入或卖出碳配额之和等于0;(4)可再生能源发电商实际售出的绿证数量不超过自身可供售出的绿证数量上限;
3.1 基础数据假设电力批发市场中共有3家发电商并且各发电商仅有一个发电机组。发电机组1和2为火力发电机组。发电机组3为可再生能源发电机组。各发电机组的相关技术经济参数可以参考文献[10]。另外,假设某时段电力批发市场逆需求函数的参数A=250,B=0.07。碳交易价格逆需求函数关系式中的参数 a 、 b 分别为40、0.05。绿证交易价格的函数关系式中的参数ζ、τ分别为300、0.02。单位时间内用户电力需求特性参数Q取2 000 MW·h。3.2 结果分析